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Abstract 

T^ We study the relation between 5D super Yang-Mills theory and the holographic de- 

scription of 6D (2, 0) superconformal theory. We start by clarifying some issues related 
to the localization oi J\f = I SYM with matter on S^. We concentrate on the case 

^S^ of a single adjoint hypermultiplet with a mass term and argue that the theory has 

a symmetry enlargement at mass M = l/(2r), where r is the S^ radius. However, 

CN in order to have a well-defined localization locus it is necessary to rotate M onto the 

imaginary axis, breaking the enlarged symmetry. Based on our prescription, the imagi- 
nary mass values are physical and we show how the localized path integral is consistent 

^^ with earlier results for 5D SYM in flat space. We then compute the free energy and 

the expectation value for a circular Wilson loop in the large N limit. The Wilson loop 

^T calculation shows a mass dependent constant rescaling between weak and strong cou- 

(y-^ pling. The Wilson loop continued back to to the enlarged symmetry point is consistent 

^~~i with a supergravity computation for an M2 brane using the standard identification of 

^ the compactification radius and the 5D coupling. If we continue back to the physical 

'k^ regime and use this value of the mass to determine the compactification radius, then 

;^ we find agreement between the SYM free energy and the corresponding supergravity 






> 



a 



calculation. We also verify numerically some of our analytic approximations. 



1 Introduction and main results 

The G-diinensional (2, 0) theories are the most mysterious in Nahm's classification of su- 
perconformal field theories [1]. They are maximally supersymmetric with no dimensionless 
parameter that can be varied between weak and strong coupling. Their fields are comprised 
of self-dual tensor multiplets and do not have a Lagrangian description. They are conjec- 
tured to be dual to M-theory on an AdSj x 5*^ background, which reduces to supergravity 
in the large N limit. One of the most striking results learned from the supergravity duals is 
an N^ dependence in their free energy and conformal anomaly [2l |3] . 

Compactification of the (2, 0) theory on a circle with radius Rq reduces it to 5-dimensional 
maximally supersymmetric Yang-Mills (MSYM) [U [S] , with the coupling identified with the 
radius by [j 

Re = ff. (1.1) 

This relation follows from the identification of the Kaluza-Klein modes in the (2, 0) theory 
with the instanton particles in the 5D MSYM [5]. 

Recently, it has been argued that the 5D MSYM can be used to define the (2, 0) theories 
[TJ [HI IS]. At first blush this seems rather strange, since 5D MSYM is not renormalizable and 
hence would appear to need extra degrees of freedom for a UV completion. In fact these 
extra degrees of freedom would seem to account for the A^^ dependence in the (2, 0) free 
energy and anomaly. In [7j it was proposed that perhaps the 5D MSYM is actually finite, 
negating the need for extra degrees of freedom. If this were true then it would be consistent 
for the 5D theory to have all degrees of freedom found in the (2, 0) theory. However, a recent 
study shows that the 5D MSYM is explicitly divergent at six loops [10] and hence requires 
a UV completion. 

Nonetheless, in [TTl [12] it was shown using localization techniques that certain 5D A/" = 1 
SYM theories on S^ with 8 supersymmetries exhibit A^^ behavior in their free energies. The 
necessary ingredient for the A^^ behavior is the presence of a single hypermultiplet in the 
adjoint representation |12], where the hypermultiplet mass sets the overall coefficient in front 
of A^^. We will call these M = 1* theories. 

In flat space a single massless hypermultiplet reproduces the field content of the 5D 
MSYM. Adding the hypermultiplet increases the number of supersymmetries from 8 to 16 
and enhances the i?-symmetry from SU{2) x f/(l) to 5*0(5). But in the Euclidean version of 
this theory on M^, the i?-symmetry must itself be rotated to the noncompact S0{4, 1) in order 
to preserve the reality conditions for the spinors. This follows from the same dimensional 
reduction argument used in the study of 4D A/" = 4 SYM [6]. 

However, unlike 4D A/" = 4 SYM, 5D MSYM is not conformal, hence there is no canonical 
way to map it from M^ to S^. In [13] it was shown how to put A/" = 1 SYM with arbitrary 
hypermultiplet content on S^ by adding additional terms to the Lagrangian. On S^, the 



"'^We normalize the Yang-Mills Lagrangian as in p], which differs from the standard normalization by a 
factor of 2, and changes the relation between Rq and gvM accordingly. 



noncompact i?-symmetry of an A/" = 1 theory is broken to S0{1, 1) x U{1). For the theory 
with a massless adjoint hypermuhiplet there is still only A/" = 1 supersymmetry. However, 
the authors in [H] observed that the global symmetry is enhanced when the hypermultiplet 
mass is M = l/(2r), where r is the S^ radius. They also provided evidence that the 
supersymmetry is increased to 16 supersymmetries, suggesting that the enhanced global 
symmetry is an i?-symmetry. A further argument is that by inserting this mass into the 
expression for the localized path integral there is a vast simplification [TT], reminiscent of 
4D J\f = 2* SYM on S"^, where the determinants in the path integral take a very simple form 
and all instanton factors cancel out at the M = A poinlQ 

Ultimately, we want to compare the free energy from 5D SYM with an analogous com- 
putation for the (2, 0) theory using its supergravity dual. It is here where the situation is 
problematic, namely because there does not exist a Euclidean version of (2, 0) theory. The 
argument for this is simple. In the Lorentzian case one can have 16 real spinors by com- 
bining an 5*0(1,5) Weyl spinor with the spinor of the 5*0(5) i?-symmetry. However, if we 
attempt a Euclidean rotation then charge conjugation maps an 50(6) spinor to the other 
spinor representation, meaning that we have to split the i?-symmetry spinors. But this is 
only possible if we split the 5*0(5) i?-symmetry into a noncompact version of SU{2) x 5f/(2). 
What we are left with is the Euclidean version of the 6D (1, 1) theory, which can be reduced 
to 5D MSYM, but with a completely different relation between the compactification radius 
and the 5D coupling. 

We can then consider the (2, 0) theory after conformally mapping it to Mx 5*^. It is natural 
to identify this 5*^ with the 5*^ for the 5D SYM. However, the M is time-like and requires a 



Euclidean rotation in order to compactify it and identify it as in (1.1). Furthermore, after 
restricting to 16 supersymmetries to preserve the 5^, the reduced /^-symmetry of the (2, 0) 
theory is 5f/(2) x f/(l), which also requires some sort of Euclidean rotation in order to 
identify it with the non-compact i?-symmetry in Euclidean 5D SYM. So it would seem that 
we would run into the same sort of problem with the spinors. 

Our approach to the problem is to just go ahead with the Euclidean rotations and see 
what we get. What we find is that we can reach agreement between the free energy of the 5D 
theory and the supergravity computation provided that we: 1) allow for a strong coupling 
renormalization of the coupling that we will explain below; 2) rotate the mass parameter M 
to i/{2r). In fact, strictly speaking it is necessary to rotate the mass onto the imaginary axis 
in order to localize the path integral. 

Localization is a powerful technique that reduces the path integral to a matrix model 
[m [TSl Uni [TT]. Using localization, the M = 1* free energy at strong coupling is found to be 

m 



9671 

where A = gyj^^N /r is the 't Hooft parameter and m = i rM. Given the scaling behavior of 
A, we see that F scales as A^^. This result can be compared to a supergravity calculation on 



^We thank V. Pestun for a discussion on this point. 



AdSj X 5''* with the AdS^ boundary chosen to be S^ x S^. In this case, using the techniques 
in [E], [TTJ UHl [19] , the regulated total action is found to be 



lAdS 



57tRq 

' 12 r 



N-^ 



;i.3) 



Using the identification in ( 1.1 ), we see that there is a mismatch of 81/80 between F and lAds 
at r/i = 0, or 4/5 if we continue (1.2) back to the symmetry enlargement point at m = i/2. 



We emphasize that we are comparing strong coupling results, so this mismatch is not like 
the case of A/" = 4 SYM at finite temperature with its famous factor of 3/4 pO] . 

While localization can only be applied in a limited number of situations, one of these is 
a supersymmetric Wilson loop which wraps an S'^-fiber of S^, where the S^ is seen as the 
Hopf vibration over a CP^. As an example, we can consider a supersymmetric Wilson loop 
that wraps the equator of the S^. At weak coupling one can show that 



(W) ~ exp 



A 

8^ 



;i-4) 



In this paper, we will show that at strong coupling the Wilson loop behaves as 

A 



{W) ~ exp (9/4 + m 



87r 



;i.5) 



In this paper we will also compute the regularized circular Wilson loop in supergravity, which 
is found by wrapping an M2 brane around an 5*^ and attaching it to a great circle on the S^ 
boundary. Here we find that 



/Tin f2nNRe 

[WjAds ~ exp 

V r 



;i.6) 



Hence, at m = i/2 we see that (1.5) and (1.6) are consistent with the identification in (1.1) 



This still leaves the mismatch in the free energy. However, if we rotate back to m = 1/2 



and use this value to match Rq to the coupling in (1.5) and (1.6), we have 



Rr 



^9ym 
327r2 



;i.7) 



Substituting this value of Rq into (6.14 we then find agreement with the supergravity com- 
putation! The mass parameter can be thought of as the expectation value of a real scalar 
field that is part of a vector multiplet. Localization reduces the path integral to a matrix 
integral over N real scalars that are also part of vector multiplets. But convergence of the 
integral requires a Euclidean rotation for all of these scalars. So perhaps it is not surprising 
that we should also rotate the mass parameter. 

The rest of the paper is organised as follows. In section [2] we review the formulation of 
J\f = 1 5D Yang-Mills theory on S^ and discuss its symmetries. In section [3] we present the 
details of the matrix model resulting from localization and explore its different limits. In 



section |4] we consider the large A^ behavior of the A/" = 1* theories. Here we calculate the free 
energy as well as the expectation value of a supersymmetric Wilson loop in the weak and 
strong coupling limits. We also generalize these results to a Z^ quiver theory. In section |5] we 
collect some numerical results about this model, showing that they agree with the analytical 
results of the previous section. In section [6] we describe the supergravity derivation of the 
free energy and Wilson loop, in particular showing that the Wilson loop is consistent with 
our identification of R^ to gy^- In section [t] we discuss some open issues. 

2 A/" = 1 5L> Yang-Mills with matter on S^ 

In this section we review the construction of A/" = 1 supersymmetric Yang-Mills theory with 
matter on the five sphere S^. This model has been constructed in [13] and we follow their 
conventions. 

Let us start with the discussion of the vector multiplet. Preserving 8 supercharges one 
may construct the A/" = 1 theory with the following Lagrangian density on 5*^ 



■'vector 2 

9ym 



Tr 



2 2 r r^ 



+z\iT"'D^\' - \j[cr, \'] - '-t''\iXj\ , (2.1) 

where we have chosen to write the radius r of S^ explicitljQ This J\f = 1 Yang-Mills theory 
is not conformally invariant and requires some guess work to construct the theory. If we 
integrate out the auxiliary field D and consider the bosonic part of the Lagrangian 



■'vector 2 

9ym 



^F^„F™" - DrnaD"'^ - 4^' + 
2 r 



2 



(2.2) 



we find an r-dependent mass-term for a. The Lagrangian density of massless scalar in 
curved space is given by 

Lscalar = D^<f)D'^<f) + |7^^^^ 0' , (2-3) 



which is invariant under the Weyl transformation of the metric Qmn -^ G^^Qmn and the scalar 
field (f) — )■ e^~^(l). Here 7^ is the scalar curvature, which for the (i-sphere is 7^ = ^ . Thus 
restricting to the case of a five dimensional sphere we get 

Lscalar = D^'pO'^'P + |^ </'' ' (2-4) 



''The unusual sign for the a kinetic term foUows the convention in [B]. We ultimately want to consider 
Euclidean Yang-Mills and yet work with physical fermions. This is accomplished by making a time-like. This 
can be seen directly for the J\f — 2 gauge theory in flat space M^ which can be reduced from lOD Yang-Mills 
on E^'^. The scalars in the 5D theory correspond to the gauge field components along the compactified 
directions. Choosing one compactified direction to be time-like so that the remaining theory is Euclidean 
results in a wrong sign kinetic term for that scalar. 



Hence, the vector multiplet scalar a is massive. 

Next we would like to couple the vector multiplet to a hypermultiplet in representation 
R. The massless hypermultiplet is conformal so there is a well-defined prescription to put it 
on the sphere. The Lagrangian density on 5*^ for the hypermultiplet coupled to the vector 
multiplet is given by the expression 

15 
Lraatter = t^'^ D^qiD'^qj - t^qja'^qj + -^e^qjqj - 2ii;Ipip - 2ipaip 

-Ae''^\iqj-iqjD'\j , (2.5) 



which contains the conformal mass term in (2.4). A more general mass term can be generated 



through the standard trick of coupling the hypermultiplet to an auxiliary vector multiplet 
and giving an expectation value to the scalar in the multiplet. This then leads to the mass 
term [13] 

Lmass = -Mh'^qjqj + -Mt'^qjqj - 2M7/>^ . (2.6) 

r 

Since a vector multiplet scalar is real, M is assumed to be real. As discussed in [121 [U] the 
localisation of the model will require the rotation of the scalar a — t- za and mass M — >■ iM, 
otherwise we will fail to get the correct localisation locus. For example, if we rotate a and 
do not rotate M then we will be unable to argue that the model is localised at g = 0. Thus 
one can construct on S*^ an TV = 1 supersymmetric Yang-Mills theory coupled to a massive 
hypermultiplet in representation R. Due to the reality conditions for the hypermultiplet, 
the representation R and its complex conjugate R enter the construction symmetrically. For 
further details on general A/" = 1 theories on S^ we refer the reader to [I3] . 

It is natural to ask if one can construct a theory on S^ with more supersymmetries. In 
fiat space the A/" = 1 theory is enhanced to A^ = 2 if there is a single massless adjoint 
hypermultiplet. But by itself there is no canonical way to map the A/" = 2 on fiat space to 



the sphere. Our best bet is to look at the scalar mass terms in (2.6) and look for a value 



for M where there is enlargement of the global symmetry. Combining the relevant terms in 



(2.4) and (2.6) we have 



4 a , 15 , ,ox rr 2i 



- i^' + (7^ - M')e'-'qjqj + -Mt'^Mj . (2.7) 



One can see that there is a special point M = ^ (or M = — ^) where the terms (2.7) become 

ED 



4 2 3 _ 1 4 _ 2 /o c 



1_ 
2r 



thus indicating an enlargement of the symmetry to 5*0(1, 2) x 5*0(2). For the point M - 
the roles of q^ and g^ are interchanged. This is the largest possible symmetry for the sphere. 
The fact that the symmetry enlargement happens for the massive hypermultiplet is because 
the vector multiplet scalar itself is massive. 
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3 Matrix model for A^ = 1 5D Yang-Mills with matter 

Using the model discussed in the previous section the perturbative partition function was 
derived in [15] for massless hypermultiplets (see also |21]). Here we discuss some subtle 
issues in the definition of the matrix model for A/" = 1 supersymmetric Yang-Mills theory. 
We follow the conventions in ^S\, spelling them out when necessary. From now on we absorb 
the radius r into the integration variable (f) = —ira. In complete analogy to 0, we must 
integrate over imaginary a, and hence real 0, in order to have a well-defined path integral 
(see |T5] for further details). 

Consider the theory with a semi-simple compact gauge group G. We have an A/" = 1 vec- 
tor multiplet and an A^ = 1 massless hypermultiplet in representation R (half-hypermultiplets 
in representations R and R). The partition function of this gauge theory on S^ is given by 
the following expression 



ZTXoM)ZrZM + 0{e 4m ) , (3.1) 



-looplrJ^l-loop 
Cartan 

where the one-loop contributions are given by the following infinite products 

a. ,1.2 



and 

^iX('^)=nn(^-(^'^'/^)+|j ■ (3-3) 

Here /3 are the roots while /i are the weights in R. In this expression everything is taken 
from ^^ except for the Tr0^ term. Let us explain its origin and its numerical coefficient. 
We define the gauge theory with matter and Chern-Simons terms 

S=^ /Tr(FA*F) + --- + -^ fTi{AAdAAdA) + --- , (3.4) 

9ym J 247r J 

where we have only written terms relevant for the normalisation. The supersymmetrization 
of the Chern-Simons term was discussed in [13]. Using the normalisations and conventions 
in [T5], in particular the relation ndndn = — 8volg where n is the contact form, we get the 



expression (3.1). 



3.1 Renormalization of the matrix model 

Ignoring for the moment the Lie algebra structure, the building block for a one-loop contri- 
bution is given by the following infinite product 

00 
V = xl[{t + a;)(i+i*+5*^) (t _ a;)(^-i*+^*') . (3.5) 

t=i 



As it stands, this infinite product is divergent. The divergent piece is extracted from the 
following expression 

logP = y^ I 3x I + convergent part . (3-6) 

Thus the relevant function can be defined by the following Weierstrass type formula 

S,(.) ^ 2.e-<''->.e«-i^n ( (l + ^)'"''^^"' (l - ff'-*''^" e*-) , (3.7) 



i=l 



where the divergent piece is subtracted from each term. Indeed the above expression is the 
definition of the triple sine function in terms of an infinite product [22]. The use of triple 
sines for 5D partition functions has been pointed out in [23], l2l] . 

Alternatively, we can introduce a cut-off in order to regularise the divergence. Cutting 
the mode expansion of the divergent part off at Uq = vrAr, the regularized 1-loop contribution 
becomes 

log(^nop(0)^i'Xp (</>)) = -^ J2(^P,^<f>)r + ^ E((^<^""))' + convergent part 

= vrAr (C2(adj) - C2{R)) Tr(02) + convergent part ,(3.8) 

where Tr(T^TB) = C2{R)Sab and Ylii^'Pi lAY = 2C2(-R)Tr(0^). We use the conventions that 

C2{R) = 1/2 for the fundamental representation of SU{N) . The linear piece disappears 
since the gauge group is semi-simple. We see that the divergent piece is proportional to 
Tr(0^) and thus can be absorbed into a redefinition of the coupling constant by 

-^ = -^--^(C2(adj)-C2(i?)) . (3.9) 

The above renormalisation of Yang-Mills coupling agrees with the known fiat space results. 
This renormalisation has been explicitly derived using supergraph techniques in [25]. Al- 
ternatively it can be derived by matching 4D and 5D prepotentials on M^ x S^ following 



jX 



Hence, the convergent part of the infinite product (3.5) can be replaced by S'3(x)e~^ 
up to irrelevant (x-independent) constants. This additional exponential factor can be ab- 
sorbed into the Tr((/)^) term by a further shift in the coupling constant, which is equivalent to 
redefining the cut-off by a finite shift. Therefore, we conclude that we can write the matrix 
model using the effective coupling g^fj and the triple sine functions 5*3. We can also see that 
in the case when C2(adj) = C2{R), the one- loop answer is automatically finite and does not 
require any regularisation. 



3.2 Massive hypermultiplet 

Assuming that infinite products are regularised (if necessary) we can rewrite our matrix 
model as follows in terms of triple sine functions 6*3(2;) 

Z = j d(t>e^' ' ' ' 'detAd[SS^))dei-R\sAict>+-\) , (3.10) 

where from now on we assume that qym is the renormalized coupling. The triple sine function 
5*3 (x) has the following symmetry properties, 

Sz{-x) = S^{x + 3) , sAx+^^= Ss (-x + ^] . (3.11) 

Assuming the standard normalization for the group generators, [Ta,Tb] = ifJ^Tc, the 
weights are switched from /x to —^ when exchanging representation R with R. Hence, we 
have the following property for the one-loop contribution of a massless hypermultiplet 

det^(^3 (# +1))= deti?(53 (-1^ +l))= ^^^r{Ss (i^ +1)) ■ (3-12) 

Therefore, the representations R and R are automatically symmetrized in the determinants, 
as is required for a hypermultiplet representation. 

Masses for hypermultiplets can be turned on easily by using the auxiliary U{1) vector 
multiplet discussed in the previous section. We simply take a G x [/(I) matrix model, 
but exclude the integratation over the U{1) direction. Thus the contribution of massive 
hypermultiplet is given by 

Z= d(l)e ^^ ' detAd[S3m) det^' i^S-i U0 + 2m+ - j j , (3.13) 



where m is related to the hypermultiplet mass parameter M in (2.6) by ?Ti = —irM. As we 
have stressed in section [2] the rotation of M to imaginary values is required by localization. 
Using the triple sine symmetry properties, we find the relation 

detR(^sJi(j) + im + -]^ =detR(^S3(t(f)-im+-]y (3.14) 

Thus the partition function with a massive hypermultiplet can be also written as 

f j^ -#^Tr(</,2)_z^Tr(^3) /^,.^A 

d(p e 'YM dei Ad[b3[^(p)j 

X det;^'/' (^83 N0 + zm + ^ n det-y^ [s^ N0 - ^m + ^ j) . (3. 15) 



3.3 Large volume limit 

Let us now study the matrix model in the large volume limit for the S^. Let us write the 
matrix model in the following form 

/ d(j) e-^ , (3.16) 

where 

^=^Tr(02) + ^Tr(03)_ J]log53((z0,/3)) + J]log53(R/x) + zm + ^) . (3.17) 

9yM ^ ^ A' 

We can restore the S^ radius dependence by the rescaling — )■ r0 and m — )■ rm. Using the 
asymptotic expansioiy for |Imz| — )■ oo and < Rez < 3 

log ^3(2) sgn(Im2;)7ri (7-2^- -z^ + z + ...j (3.18) 

we obtain the following behaviour 

^^ = ^Tr(0^) + ^Tr(03) + ^ ($^ |(0, /3)|^ - $^ |(0,/i) + mf) + 0{r-') . (3.19) 

Modulo a constant which was absorbed into the definition of the coupling in [271 128] , this 
expression matches the exact quantum prepotential in the fiat-space limit |2H] ■ The normal- 
isation in front of the quadratic term can be fixed either by a direct one-loop calculation in 
fiat space |25] or by matching the 5D superpotential with the 4D superpotential as in [26] . 



Notice that M in (2.6) must be rotated to the imaginary axis in order for m to match the 



mass parameter of the fiat space physical theory [28] . 

3.4 Well-defined matrix model 

We next ask under what conditions the matrix model is well-defined, i.e. the matrix integral 
converges. We answer this by going to large (p and finding if J-" is a convex positive function 
in this limit. Taking the limit, we get the following asymptotic expression 

9ym 



^=^Tr(0^) + ^Tr(0^) + ^(5^|(0,/3)|^-$^|(0,/.)|= 

9y1vI ^ ^ 

-|m5^sgn((0,/.))((0,/z))^-7r5^|(0,/3)|-|(m2 + iW|(0,/.)| + --- , (3.20) 

^In this asymptotic expansion the dots stand for the constant term. This expansion is typically derived 
for the region < Rez < 3 using the integral representation for 5*3(2:). The case Rez = is derived separately 
via the explicit expression of 6*3 in terms of polylogs. 

10 



where ■ ■ ■ stands for terms suppressed for large 0. Analyzing the convexity of the above 
function, it is clear that the cubic terms dominate. Hence, the analysis is identical to the 
one presented by [28] and so the same conditions apply. In some special cases the cubic terms 
can cancel each other. For example, this happens in the case of single adjoint hypermutiplet 
[12] or in the case of USp{2N) theory with specific matter content [29] . 



3.5 Decoupling of a massive hypermultiplet 

We now consider the behavior of the matrix model as we send the hypermultiplet mass to 



infinity. For large \m\ the leading terms in (|3.17|) are 



^vr^r 



Tr(02) + !5^Tr(03) - J]logS3((^0,/3)) -sgn(m)^ J] (l((0,^))=' + m((0,/i)) 



9ym 3 2 \ 3 



M 



The two last terms can be absorbed by a redefinition of k and Qym- 
To see this explicitly, we note that 

TiiTATBTc + TaTcTb) = CsiR)dABC , (3.21) 

where we used the following relation 

5^((0,/x)f = C3(i?)Tr(03). (3.22) 

The coefficient C3 satisfies C^lR) = —CslR), hence it is zero for real representations. For 
the lower complex representations in SU{N) it is 1 for the fundamental, iV — 4 for the 



antisymmetric, and iV + 4 for the symmetric representations. Hence, from (3.21) and (3.22) 
we get 

fce// = fc-sgn(m)^l^, (3.23) 

reproducing the result in [271 EH]- Similar analysis of the quadratic terms gives the formula 

4 Af =r 5D super Yang-Mills 

From now on we concentrate on a single adjoint hypermultiplet with mass parameter m, 
which we refer to as A/" = 1* super Yang-Mills. We also assume that the gauge group is 

SU{N). 
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4.1 The free energy 



In order to analyse the matrix model (3.15) we use an alternative representation of the triple 

(4.1) 



sine 5*3(2;) as defined in (3.7). Namely we can rewrite Ss{z) as 



S3{z) = 2e-^'(-2) sin(7rz) e^^'^ ei'(") , 
where l{z) is the function introduced by Jafferis [30] 

l{z 

and f{z) is the function introduced in 



\og{l-e'-^^) + '-{nz' + hi,{e'-^^^ 



TX 



ITX 

12 



ZTTZ 



IZ, 
7T 



1 



C(3) 



+ z^ log (1 - e-2-*^) + -Li2 {e-^""") + — ^Lig (e-^--) 



(4.2) 



(4.3) 



Using this representation the matrix model (3.15) can be rewritten as follows 

]^(sin(7r(/3,^0))e-3'(i-^— </^'*<^))-3'(i-^™+</5'*<^)) 



Z = / [ci0] e 

Carian 






X g|/{(/3,i</.»-i/(|-im-(/3,i0»-i/{|-im+(/3,i</,» _^ 



(4.4) 



where ■ • ■ stands for the contribution from instantons. Our main focus is on Yang-Mils, so 
we have dropped the Chern-Simons term. Using the t' Hooft coupling constant 



A 



91mN 



and taking the large A^ limit at fixed A, all instanton contributions are suppressed and the 



partition function in (4.4) reduces to the matrix integral. Specialising to SU{N) we rewrite 



the partition function (4.4) in terms of the eigenvalues (pi and end up with the following 
matrix model 



Z 



N 

n 



;exp 



A 






log [sinh(7r(0i - 0^))] 



--/ i--im + i{q 



HI 



im - i{^i - 0j) ) + -/(i(0i - 0j)) 



-fi im — i( 



(4.5) 



In order to proceed further let us review the relevant properties of the /- and /-functions. 
These are: 

1. l{z) is an odd function and f{z) is an even function 

l{z) = -l{-z), f{z) = f{-z) 
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2. The derivatives of the functions are given by the following simple expressions 

dl{z) 



; = TTZ COtfTTZ) 

dz 



dz 



—nzcoti-KZ] 



(4.6) 



3. The asymptotic behavior of the functions are given by 



lim Re/ I — \- ix 

|a:|— s>oo \ 2 



■-\x\^ + -\x\; lim^^^lmf [^±ix) =±^^^ , 



-'^rx^ 



lim Re/ I — \- ix 

|x|^oo V 2 



IT 

"2 



x\ ; lima;^oo Im/ (^ ± ix) = =Ff x 



2 . 



(4.7) 



lim Ref (ix) = — — \x\" ; 

|a;|— ^oo o 



Im/ (ix) = . 



In the large N limit the partition function is dominated by the saddle point. Using the 



derivatives in (4.6) the (pi satisfy 



X 






(2 - {<j)i - <j)^y) coth(7r(</., - </.,)) 

— (f)j — mY J tanh(7r(0j — (pj — m)) 



+ 2(4 + 

1 /I 
+ 2 4 + 



— (j)j + mf j tanh(7r(0j — (pj + m)) 



For weak coupling where A <^ 1, (|4.8|) reduces to 



IQti^N 



whose solution has the usual Wigner distribution 



P(0) 



1 dn 



N d(p TTcpl ^ ^° 



^2 _ A2 



47^3 



where 



p{<P)d<P = 1 . 



In this limit the free energy has the typical weak coupling form 

F = -logZ^ -A^^log^A. 



(4i 



(4.9) 



(4.10) 



(4.11) 



(4.12) 



13 



In the opposite limit where A ^ 1, we can assume that 10, — 0^ | ^ 1 if we further assume 
that |m| ^ A. In this case the saddle point equation simplifies to 






9 

TT I — \~ m" 
4 



Y^ sign( 



(4.13) 



j¥=i 



Assuming that the eigenvalues 0j are ordered, we get the solution 

(9 + 4m2) A 



which corresponds to an eigenvalue density 

327r2 
p(0) = ^— \d\ < 



i2z~N), 



(4.14) 



(9 + 4m2) A 




(9 + 4m2) A 
64^^2 



> 



(4.15) 



Using the asymptotic expressions in (4.7), the matrix model (4.4) at strong coupling can 
be simplified to 



n 



he 



^E0?+f(|+-^)EEI<^.-</'.l 



(4.16) 



Substituting the saddle point solution (4.14) back into (4.16), we find the free energy, 

F = -lofyZ^- ''i':'' (- + m^l , (4.17) 



where we used the approximations 



N N 

^(2^ - iV)2 ^ -iV^ ^^|^-j-|^-iv3. (4.18) 

j^i i=l 



i=l 



Hence, we see that going from weak to strong coupling the free energy crosses over from A^^ 
to N^ behavior. At the point m = ^ the free energy is 

F = ^^^^ — . 4.19 

3847rr ^ ^ 

It is interesting to see how the system evolves as m approaches a large value. In this 
limit we expect the system to reduce to A/" = 1 SYM with no hypermultiplets. Indeed, if we 
assume that m :» |0j — 0j| for all i and j, the saddle point equation reduces to 



leTT^A^ 






(2 - (0, - 0^)2) coth(7r(0i - 0j)) + 2m (0i - 0^) 



27imN(pi + TT ^ (2 - (0i - cf>jf) coth(7r(0i - 0^-)) 



(4.20) 
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where we used that J2i^i — 0- We can then reexpress (4.20) as 



TT J] (2 - (0, - 0,)2) coth(7r(0, - 0,)) , (4.21) 



^eff ... 

which is the saddle point equation for the A/" = 1 system with no hypermultiplets and an 
effective 't Hooft couphng 

11m , , 

(4.22) 



^eff A 



ill 



2 



Note that this is consistent with our consideration in section |3] and the exphcit one-loop fiat 
result where m is regarded as UV-regulator. 

4.2 Supersymmetric Wilson loops 

Supersymmetric Wilson loops can be evaluated using localization since they preserve some 
of the supersymmetries. In 5D flat space such loops were first considered in ^T\ . Supersym- 
metric Wilson loops were also considered in [521 133] • 

On 5*^, the supersymmetric loop must go along an S^ fiber when the S^ is viewed as 
an S^-fibration over CP^. Thus proceeding in analogy with [6] the expectation value of 
supersymmetric Wilson loop is given by the following matrix model expectation value 

(1^) = i-(Tre2'^'^0 . (4.23) 

Taking into account the consideration from previous subsections we first observe that in the 
large A^ limit the Y2 e'^'^'^' term has a negligible back-reaction on the position of the saddle 



point in (4.14). Thus at strong coupling A we have to evaluate the integral 



(W) -^ - lld<p,J2e'-'^^e ' ^>^- . (4.24) 

•^ i i 

Thus at large A^ the Wilson loop expectation value is well-approximated by the following 
integral 

{W) = j d(Pp{(P)e^'"^ (4.25) 



where p(0) is the density of eigenvalues. 

For weak coupling we can approximate the integral as 

{W) ^ /c/0p(0)(l + 271^02) = 1 + A ~ exp ( -^ ) . (4.26) 
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At strong coupling, using the eigenvalue distribution in (4.15), we find 



where we have omitted the prefactor, as only the exponential term is important for us. 
Comparing the weak and strong coupling results, we see that at large A the effect on the 
Wilson loop is to rescale A by the factor (| + m^). 

4.3 Quivers 

The above results can be straightforwardly generalized to a Z^ quiver of the A^ = 1* the- 
ory [12]. Here we have an SU{N/kY gauge group and equal mass hypermultiplets in the 
bifundamental representations, {N/k, N/k, 1, ... 1), (1, N/k, N/k, 1, . . . ), etc.. The eigenval- 
ues of (4.8 ) are split into k groups V^j , where r = 1, . . . , fc, -i = 1, . . . N/k, and the resulting 



-# = 7r[5^(2-(^f)-^f))2)coth(vr(^f)-^f))) 



equations of motion are 

levr^A^ 

•pr = T^[ 

+ (E [i (l + (#-^r^-^)0 tanh(vr(#-^f ^)-m)) 
+ i (I + (#-V^r^-m)^) tanh(vr(#-^f-^)-m))]^ 

+ (m -> -m) I . (4.28) 

(r) (s) 

These have a solution where ip^ = Wi or all r and s, hence taking the same limits as before 
we find that the eigenvalues sit at 



(9 + 4m2) A 

Wi ' = 

Thus, the free energy is 



'^'-m^^^'-'^i'^y (^-2^) 



F--^'M^iUrA'-%^iUr^\\ (4.30) 



967rrA;3 y4 J 967rrk'^ y4 

For supersymmetric Wilson loops, the weak coupling behavior parallels the A/" = 1* case 
with A replaced by X/k, hence 

(,y)«exp(^). (4.31) 
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Figure 1: Density of eigenvalues p(0) for ?7i = 0, A^ = 200, /3 
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N = 160, 13 
solutions. 



9ym 



2' 



2 (right). The dashed blue lines are the analytical strong coupling 



At strong coupling the leading behavior of the Wilson loop is determined by the top eigen- 
value, hence we have 



(W) ~ exp 



A /9 



+ m" 



871 k V4 
Therefore, we find the same large A rescaling of the Wilson loop as in the A/" = 1* case 



(4.32) 



5 Numerical study oi Af = 1 bD Yang-Mills 

In this section we give numerical evidence for the analytical approximations used in the 
previous section. 



We have been unable to find an exact solution for the saddle-point equation in (4.8) 



However, we can look for numerical solutions using an idea similar to the one used in [M] . 
Instead of solving a system of A^ algebraic equations (4.8) of the form — ^ = 0, where J-" 



is defined in (3.21), we introduce a time dependence for the matrix model eigenvalues </'i(t) 



and solve the "heat" equation 



dJ^ 



T- 



dt 



(5.1) 



At large time-scales where t — )■ oo, with an appropriate choice of r the solution of (5.1) 
relaxes and approaches the solution of the saddle point equations — |x" = 0. 

Following this approach we eventually reach the density of eigenvalues p(0) = ^ shown 
in Fig. [1} Here we show the densities for two different values of m. For comparison, we have 
superimposed these over the corresponding analytical strong coupling result from (4.15). As 



one can see the analytical and numerical solutions coincide at large A^, except for a small 
region near the boundaries of the distribution. In the case of pure A/" = 1 SYM theory with 
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Figure 2: Density of states for pure J\f = 1 SYM with (3 = 2. 
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Figure 3: (^/e/ijA^-dependence of free energies for different values of hypermultiplet mass m: 
m = orange., ra = \ purple, m = 1 blue; (right) Free energy A^-dependence for pure 
TV = 1 SYM. {13 = 2 for all plots.) 

no liypermultiplets the density of states develops a peak at each end as well as a shallow 
minimum in the center. This is shown in Fig. |2} 

Using these distributions of eigenvalues we can find the free energy for the matrix model. 
Results for the energies with different values of m and for TV = 1 are shown on FigjS) The 
dashed lines correspond to the respective strong coupling results from (4.17). As the graphs 



show, the difference between the analytic and numerical results grows as m decreases but 
converges as A^, and hence A is increased. The discrepancy we observe is due to the effect of 
sub leading terms in A^ and A. 

Having the data for the different m we can fit the A^-dependencies with a simple polyno- 
mial function F = a^N'^ + asN^ + 02 A^^ + aiN + qq + •••• This fit eventually gives the leading 
A^^ term for general m. The coefficient 03 can be fit to a polynomial function of m, 



as = ci + C2m + c^m^ + C4m^ + 



(5.2) 



Eventually we end up with coefficients that are very close to the strong coupling result in 
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Figure 4: m-dependence of A^^ coefficient in the matrix model free energy. 



(4.17). In Fig. we plot the ?7i-dependence of the N"^ coefficient of the free energy. The 
dashed line shows the function in (4.17). The plot shows precise agreement between the 
numerical and analytical results. 

For pure A^ = 1 SYM the fit gives the leading term behaving as A^^ instead of A^^ as in 
the A/" = 1* theories. 



6 Supergravity comparisons 

In this section we compare our strong coupling results for 5D SYM with the analogous 
computations in supergravity. To start we review the supergravity computation of the free 
energy given in [12J. We consider supergravity on AdS-j x 5^ where the AdS-j boundary is 
S^ ^S^. The radii of AdS^ and S'' are I and £/2 respectively, where i = 2lpi{-nNf/^. We 
write the AdS-j metric as 



ds"" 



' (cosh^ p dr"^ + dp^ + sinh^ p dVt 



(6.1) 



where dVt\ is the unit 5-sphere metric. The Euclidean time direction is compactified and has 
the identification r ^ r + 27ri?6/r, while R^ and r are the boundary radii of S*^ and S^. 

Under the AdS/CFT correspondence, the supergravity classical action equals the free 
energy of the boundary field theory. The action needs to be regulated by adding counterterms 
[161 IT71 ITSl 119] . There can be scheme dependence in the regulation \\3\, but we will follow 
the minimal subtraction prescription, which is the normal procedure when regulating the 
action. The full action then has the form 



lAdS — -^bulk + I'. 



AdS — -fbulk -r -(surface "T Jet : 



where 



'bulk 



167rG 



l^Yo\{S') j £x^g{R-2K) 



(6.2) 



(6.3) 
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is the action in the bulk, Jgurf is the surface contribution and I^ contains counterterms written 
only in terms of the boundary metric and which cancel off divergences in /bulk- Newton's 
constant is given by Gn = IGii'^ipi [35] . The equations of motion lead to 

12 
R-2A = --, (6.4) 

which when substituted back into the action gives 



^buik = -^^^^ [-^) -^^%-l2n J^ coshpsinh>dp = ^iV^sinhVo. 

(6.5) 

In the limit that po — ^ 00 the integral is divergent and corresponds to a UV divergence 
for the boundary theory. In terms of an e expansion of the boundary theory, we make the 
identification e = e~'^°, which then gives 



The surface term contributes to the divergent pieces, but not the finite part of (6.5), while 
the effect of the counterterm with minimal subtraction is to cancel off the divergent pieces. 
Hence, we find [17] 

/... ^ -'-^n:' . (6.7) 

We next consider the supergravity calculation [36] for the Wilson loop which is related 
to the extremized world- volume of the membrane 

(iy)~e-^''^/'^^, (6.8) 

where T^'^^ = .^Jap is the tension of the M2 brane. The M2 brane is chosen to wrap 
the Euclidean time direction and the equator of S^. The third direction falls in from the 



boundary into the bulk. Hence, the M2 brane volume is given by 

27rflg 

r 27r po 

dV = 1^ / dr d(j) dpsmh{p) cosh(p) (6.9) 



000 



Using the same UV cutoff as in (6.5) we find 



nNRe fl 



T(^) dV= ^ --2 + e . (6.10) 



As in the case of the action, the integral needs to be regulated. Using minimal subtraction 
again, the result for the regulated Wilson loop is 

(W^)~exp(?^) (6.11) 
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Comparing this to the strongly coupled result in (4.27), we find that we should set 

9ym ( 9 



Re, 



m 



167r2 ^4 

in order for them to agree. At the m = 1/2 point the relation is 

9ym 5 



(6.12) 



Rr 



167r2 2 



(6.13) 



Using this relation we see that (6.7) agrees with (4.19). 



We can also compare supergravity results for the quiver. The effect of the quiver on the 
supergravity computation is to replace the S^ with S^/Zk. This reduces the volume factor 
of the S'^ by a factor of k, hence 



12r k 



(6.14) 



Comparing the supergravity Wilson loop in (6.11) to the quiver Wilson loop in (4.32), we 
see that 



R^ = lyM^(l+rn' 



IQn^k \A 



(6.15) 



Again choosing m = 1/2 we have agreement with the free energy in (4.30) 



7 Discussion 

In this paper we have studied the matrix model which corresponds to the full perturbative 
partition function for A/" = 1 Yang-Mills theory with matter on the five-sphere. We have 
shown how the localized matrix model relates to the fiat-space results on the Coulomb 
branch in [271 EH] , [26] and explicit one-loop calculations in fiat space [25j . This indicates 
that the prescription we follow, which is analogous to the 4D prescription in [6], is consistent. 
Moreover we have shown that one can have agreement between the free energy of the A/" = 1* 
theory with a specific value of the hypermultiplet mass and the regularized supergravity 
action if one incorporates a rescaling of the physical coupling as well as a Euclidean rotation 
of the mass-parameter. The Euclidean rotation of the mass-parameter is needed to localize 
the partition function. 

An important open problem is to better understand the Euclidean version of the 6D 
(2, 0) theory, or if it does not exist, an appropriate substitute. This includes its symmetries 
and its reduction to the Euclidean 5D model, which is required for localization. This may 
provide a stronger argument for the calculations presented in this paper. 
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